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Qh, Abstract We investigate analytically and numerically a nonlinear modification of the 

2 ! magnetospheric plasma density under the action of the ponderomotive force induced by 

I ULF traveling waves, using the nonlinear stationary force balance equation. This equation 
is applied to both the dipole and dayside magnetosphere having one and two minima of the 
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G> ■ geomagnetic field near the magnetospheric boundary. The separate and joint actions of the 
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ponderomotive, centrifugal, and gravitational forces on the density distribution are shown. 
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1 Introduction 

The ultra-low frequency (ULF) electromagnetic pulsations ranging from tenths of Hz 
up to a few Hz are a common phenomenon in the Earth's magnetosphere. These waves 
are electromagnetic in character, i.e., have the electric and magnetic components. This 
is confirmed by both ground measurements (e.g. Troitskaya 1961; Tepley 1961; Heacock 
and Akasofu 1973; Bortnik et al. 2008) and in-situ observations by satellites (e.g. Young 
et al. 1981; Erlandson et al. 1990; Halford et al. 2010; Kim et al. 2010). The ULF 
pulsations are duct along geomagnetic field lines in the direction of which the parameters 
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of magnetosphere are nonuniform. This inhomogeneity must be taken into account when 
considering the influence of the wave action on the magnetospheric plasma to obtain a more 
realistic picture of physical processes. 

The theory of electromagnetic ion cyclotron waves in the ULF range has two aspects. 
The first one is the problem of sources of these waves related to the physics of the 
magnetosphere. This aspect is important not only for explaining the origin of these waves 
and mechanisms by which they are generated, but also for processes being responsible for 
the acceleration, diffusion and precipitation of magnetospheric particles. These problems 
have got a great deal of attention in past years (see e.g. Kennel and Petschek 1966; Cornwall 
1966; Troitskaya and Guglielmi 1967; Feygin and Yakimenko,1971; Gendrin et al. 1971; 
Tverskoy 1971; Nekrasov 1987; Gughelmi and Pokhotelov 1996). Another important aspect 
of the theory is the investigation of the influence of ULF pulsations on the background space 
plasma. Ponderomotive forces induced by these waves can significantly contribute to the 
plasma balance in the Earth's magnetosphere (Allan 1992; Gughelmi et al. 1993; Gughelmi 
and Pokhotelov 1994; Witt et al. 1995; Allan and Manuel 1996; Pokhotelov et al. 1996; 
Feygin et al. 1998; Nekrasov and Feygin 2005). In the papers by Gughelmi et al. (1993, 
1995), it was demonstrated for the dipole magnetosphere that a pronounced maximum of 
plasma density is formed in the vicinity of the minimum magnetic field force along the wave 
trajectory, if the amplitude of electromagnetic ion cyclotron waves exceeds a certain critical 
value. Near the dayside boundary of the magnetosphere, the geometry of geomagnetic 
field lines changes from the structure with only one magnetic field force minimum to the 
structure with two minima (so-called magnetic holes) along the field line (Antonova and 
Shabansky 1968). Pokhotelov et al. (1996) have shown that the ponderomotive force can 
lead to plasma density accumulation in the vicinity of these magnetic holes. 

In some papers cited above (e.g. Guglielmi et al. 1993, 1995), it has been assumed that 
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the equation of the force balance along the magnetic field lines contains the total thermal 
pressure together with the nonlinear ponderomotive force derived for small perturbations. 
However, it is not correct from the point of view of the nonlinear theory, in which all terms 
must have the same order of magnitude in one equation (for the correct approach see e.g. 
Allan and Manuel 1996; Nekrasov and Feygin 2005, 2011). Therefore, conclusions made in 
the corresponding papers are doubtful. 

In this paper, we consider analytically and numerically the nonlinear modification of 
the magnetospheric plasma density, using the nonlinear stationary force balance equation. 
We apply our investigation to the dipole and dayside magnetosphere having one and two 
minima of the geomagnetic field near the magnetospheric boundary, respectively. 

The paper is organized as follows. In Sect. 2, we give the expression for the 
ponderomotive force of low-frequency waves traveling along the nonuniform geomagnetic 
field. The geomagnetic field at the dayside of the magnetosphere is discussed in Sect. 3. 
The stationary force balance equation is considered in Sect. 4. In Sect. 5, numerical 
calculations of balance equation in the curved geomagnetic field are given. Conclusive 
remarks are summarized in Sect. 6. 



2 Ponderomotive force 

In this paper, we use an expression for the ponderomotive force Fp^^ along the 
background nonuniform magnetic field induced by the circularly-polarized electromagnetic 
WKB- waves traveling along field lines (see e.g. Nekrasov and Feygin 2005). In the case 
N"^ » 1, this expression has the form 
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where 

N'^^^ —r^ r (2) 

is the refractive index. Here, ojpi — {Anne^ /rrii)^^^ and oJi — eB/rriiC are the ion plasma and 
cyclotron frequencies, respectively, p — rriin is the equilibrium plasma mass density, n is the 
local number density, B is the local background magnetic field, Ei is the wave amphtude, 
uj (> 0) is the wave frequency, k„ is the wave number, e and rrii are the charge and mass of 
the ions, a denotes the left-(+l) and right-(— 1) polarization, and c is the speed of hght in 
vacuum. 

In the WKB approximation, we have Ei ~ A?""^/^ or EfN — EIqNq — const, where 
the subscript here and below relates to the corresponding values at the magnetic equator. 
We will express the ponderomotive force through the amplitude of the wave magnetic field 
BiQ at the magnetic equator because satellite measurements give usually the magnitude of 
magnetic field perturbations. Prom Faraday's equation for the traveling wave, it is followed 
that Bi = NEi. Thus, we obtain from these relations and (2) 



where vq — aui/uiiQ. Substituting (3) in (1), we obtain 
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3 Curved dayside geomagnetic field, field line equation and longitudinal 

gradient V,, 



The magnetic data collected by satellites and theoretical models show a complex 
structure of the geomagnetic field at the midday boundary of the Earth's magnetosphere. 
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This structure is characterized by a smooth transition from "V" to a "W" magnetic field 
shape with two minima in the magnetic field along field lines. There exist several analytical 
(Hones 1963; Mead 1964; Antonova and Shabansky 1968) and quantitative models (Mead 
and Fairfield 1975; Tsyganenko 1987) of the magnetic field in the Earth's magnetosphere 
based on sateUite observations. The model by Antonova and Shabansky (1968) is quite 
convenient for an analysis of geophysical phenomena in the dayside magnetosphere and, in 
addition, is in a reasonable agreement with the magnetometer data provided by the HEOS 
1, 2 sateUites (Antonova et al. 1983). In this model, the structure of the geomagnetic field is 
determined by the existence of two dipoles: an original dipole having the Earth's magnetic 
moment M and an additional one which imitates the distortion of the magnetic field due 
to the solar wind pressure. The latter dipole has the magnetic moment kM, where A; is a 
constant parameter, and is shifted a distance a (measured in the units of the Earth's radius 
Re) at the dayside along the Earth-Sun line from the position of the original dipole. 

In the spherical coordinate system, magnetic field components for the two-dipole model 
in the meridional noon-midnight plane have the form (Antonova and Shabansky 1968) 

= (5) 

= ^Vl^p^ (6) 

where x — simp, (p is the geomagnetic latitude, r is measured in units Re, and Be — 0.311 
G is the equatorial magnetic field at the Earth's surface. Coefficients a and /3 arc the 
following (Antonova and Shabansky 1968): 

a = 1 ^ (7) 

2 {a^ + - 2ar^/l^y 

kr^ [Vl - (a^ + r^) - ar {2 + x'^)] 

(3 = 1 -\ ^/ TT- (8) 

(a2 + r2 - 2arVT^^ ' 



-6- 



The value of the magnetic field B in an arbitrary point of the field fine can be defined from 

(5) and (6) as 

[ix'a' + (1 - x') P'] . (9) 

An equation for the field line is determined by dr/rd(p — B^/B^ or 

dr 2xr a 



dx 1 — x'^ (3 
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The intensity of the magnetic field force along the near boundary field fines described 
by (9) can have two minima located symmetrically relative to the magnetic equator 
(Antonova and Shabansky 1968). When a tends to infinity, the values a and /3 (see 7 
and 8) approach to 1. In this case, we have a transition to the one-dipole approximation. 
Thus, all the expressions obtained in this paper can be transformed to those corresponding 
to the one-dipole model. Since a purpose of our paper is a qualitative description of the 
ponderomotive force action induced by geomagnetic pulsations in the dayside magnetosphere 
along field lines at different distances from the Earth (up to the magnetopause) , we will not 
focus on details of the two-dipole model. 

The operator V,, is defined by relation V,, = b • V, where b is the unit vector along the 
magnetic field. Using (5), (6), (9) and (10), we find 

v„ = ^R-E'v-'^'i^, (11) 

where 

In order to investigate the ponderomotive force action induced by geomagnetic 
pulsations in the dayside magnetosphere along field lines and carry out analytical and 
numerical calculations at different distances from the Earth (up to the magnetopause), we 
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use the magnetic field model by Antonova and Shabansky (1968) described above. For the 
equilibrium plasma mass density, we take the power law form to describe the longitudinal 
field line distribution, p ~ . For the large distances from the Earth's surface which 
we consider below, the best choice for 7 to be appropriate to experimental data is 7 = 1 
(Denton et al. 2006). Thus, we take approximately 

p{x)^p^{l-x^Y\ (13) 

This formula can be applied up to </? ±50 — 60° (Denton et al. 2006). For simphcity, we 
do not take into account in (13) the non-dipole form of the magnetic field (see 10) because 
of an uncertainty of plasma mass distribution in this region. 

4 Stationary force balance equation 

From equations of motion for the ions and electrons in the second approximation on 
the wave amplitude averaged over fast oscillations, we obtain the force balance equation in 
the stationary state along the magnetic field hne (e.g. Nekrasov and Feygin 2005). Taking 
into account the gravitational and centrifugal forces (e.g. Lemaire 1989; Persoon et al. 
2009), we have 

V„P2 - {g» + nn,,fl^REr cos^ ip) p2 = Fp„. (14) 

Here p2 = P2C^ and p2 are nonlinear stationary perturbations of pressure and mass 
density, respectively, Cg — {2T/mi)^^^ is the sound speed, T is the temperature, 
g'li = g • b = —QEBr/r'^B is the longitudinal gravitational acceleration [qe — 9.8 m sec~^), 
Q is the Earth's rotation frequency, and nf^,, = nf7 ■ b, where no is the unit vector along the 
centrifugal force. Calculations show that na,, = (1 + {3 /2a) {B^/ B) cosLp. Substituting (4) 
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into (14) and using (5), (11), and (13), we find 
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dx 



(1_ ^2)1/2 

This differential equation determines a nonlinear redistribution of the plasma density due 
to the action of the ponderomotive force. 



Equation (15) can be rewritten in the form 



^^ = AA+A,iA, + A, + A,) 
dx Po Po 
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Coefficients Aj, i = 1, 2, 3, 4, 5, are the following: 
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where we have used (9) and (12) in the first term on the right-hand side of (16). The 
acceleration geff in the coefficient Ai is equal to 

geff ^gE + gcf = - (l + P/2a) Q^Ret^ cos^ if, (18) 

where gE and gcf describe an influence of the gravitational and centrifugal forces, 
respectively. 

If we take for estimation a ~ /3 ~ 1 in (18) and substitute numerical values of 
parameters gE, fi, and Re, we obtain 

gE 3 X 10^ 



gcf cos^ (fi ' 



(19) 
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where r fa rocos^f/? (see 10). Thus, for rocos^c/? < 6.7, one can neglect the contribution of 
the centrifugal force. For example, the last condition for tq = 8 is satisfied at \(p\ > 20°. 
The same for tq = 10 is at \ip\ > 30°. In the region of the magnetic equator, the centrifugal 
force is dominant and leads to a local peak in mass density. An estimation for ro — 8 
corresponds to the results by Denton et al. (2006) (see their Figs. 9 and 10). 

The coefficient Ai given by (17) will be the same in the equation for the background 
equilibrium. Its estimation for gravity, for example, at ro = 8 is the following: 

~ 8 X 10~^a; (1 — a;^)~^ (we have taken = 10^ m^ sec^^ for T = 10 eV). Thus, one 
could wait an essential contribution of the gravitational force to the equilibrium mass 
density distribution at 1 — < 4 x 10~^. However, this result does not coincide with (13). 
This could be a consequence of additional factors except for the gravitational force which 
lead to the power law distribution of mass density. 

We see from (16) that the effective gravity (?e// and ponderomotive force ~ A2 influence 
on the nonlinear density redistribution. These two forces can be of the same order of 
magnitude or one force can dominate in the dependence on the wave amplitude. We note 
also that the nonlinear disturbance of density increases when i/qBo ^ B or ou ^ uji. This 
case can occur in magnetic holes. 

5 Numerical results 

For numerical calculations, we have chosen a = 33 and k — 13. Such parameters 
correspond to the dayside boundary of the magnetosphere at the distance IORe and the 
region of experimental data obtained by HEOS 1, 2 satellites (Antonova et al. 1983). 
Solution of Eq. (16) has been carried out by means of the Runge-Kutta method. In all 
the cases considered, we assumed that the square of the sound velocity = 10^ m^ sec~^. 



- 10 - 



As it has been mentioned above, this value corresponds to T = 10 eV. The plasma mass 
density at the equator po was chosen equal to 1.67 x 10~^° kg m~^ for all L (L is Mcllwain's 
parameter) near the midday boundary of the Earth's magnetosphere. These values are 
based on experimental data from measurements of the plasma density from the plasmapause 
up to the magnetosphere boundary. In these areas, the plasma density depends weakly 
on L (Chappel 1974; Carpenter and Anderson 1992). Since the purpose of this paper is a 
qualitative study of the ponderomotive force effect of low-frequency perturbations on the 
redistribution of the background plasma density, we do not focus on details of the weak 
dependence of po on L. 

The dependence of the Earth's magnetic field B on the geomagnetic latitude (p in the 
meridional plane of the dayside magnetosphere for different L for the model by Antonova 
and Shabansky (1968) is given in Fig. 1. This figure shows a smooth transition from the 
"V" to "W" magnetic field shape with two minima in the magnetic field along the field 
lines. Fig. 2 demonstrates the separate and joint infiuence of the gravitational, centrifugal, 
and ponderomotive forces on the relative latitudinal density perturbation 6 = P2/P0 for 
L — 6 and k = 13 (see 16 and 17). The same is shown in Fig. 3 at L = 10. We see that 
the role of gravity decreases and the mass density is peaked due to ponderomotive and 
centrifugal forces. Such peaking is observed in the magnetosphere of the Earth (Denton et 
al. 2006). Figures 4 and 5 represent the difference in S distribution for the dipole {k = 0) 
and non-dipole {k — 13) geomagnetic field at L = 8 and L = 10, respectively. We see that 
the essential difference is due to the presence of geomagnetic holes. We note that small 
variations in the density distribution given in Figs. 2-5 are because of the small value of the 
coefficient Ai in (16). In order to show different dependence's, we have taken a small value 
for Bio. 

At the magnetic equator, the projection of the ponderomotive, gravitational, and 
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centrifugal forces on the geomagnetic field line is equal to zero. The longitudinal centrifugal 
force /c/ii is expressed via the longitudinal gravitational force fg„ at cos 93 ~ 1 and for 
To = 8 -j- 10 in the form 

/c/ii ~ 1.7^3.3/p„ 

(sec 19). Thus roughly, these forces are of the same order of magnitude and small (see 
discussion below 19). Their influence is displayed when the ponderomotive force is also 
small (Figs. 2-5). At the wave amplitudes Bio larger than 10"^ G, the gravitational and 
centrifugal forces are not important. 

6 Conclusion 

In this paper, we have performed an analytical and numerical study of the influence of 
the geomagnetic field structure near the dayside magnetosphcric boundary on the plasma 
density redistribution due to the ponderomotive force action induced by ULF perturbations. 
We have shown that the ponderomotive force causes the nonlinear stationary plasma density 
redistribution, which is much smaller than the background plasma density. We note that 
the ponderomotive force is a nonlinear effect, which is found by using linear perturbations. 
The assumption that the ponderomotive force has the same order of magnitude as the 
background pressure and gravity is incorrect. However, such an assumption has been 
adopted in some papers (Guglielmi et al. 1993, 1995; Guglielmi and Pokhotclov 1994; 
Pokhotelov et al. 1996; Feygin et al. 1998). Ponderomotive force (1) is derived by the 
perturbation method. Therefore, its contribution can only disturb an equilibrium state. It 
is followed from (15) that P2/P0 ~ BIq/AttpocI on the order of magnitude (see also Nekrasov 
and Feygin 2005). For parameters given above and for Biq = 5 x 10~^ G, the ratio P2/P0 is 
of the order of 10~^. 

We have investigated the plasma density redistribution for the model of the geomagnetic 
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field B by Antonova and Shabansky (1968) (Fig. 1). This choice is justified by its simphcity. 
The use of a more complex geometry for B can be useful for an analysis of certain particular 
situations (for example, specific satellite measurements etc). The results, however, will be 
qualitatively the same. We have shown that the two-dipole structure of the Earth's magnetic 
field (geomagnetic holes) leads to an essential modification in the plasma distribution (Figs. 
4 and 5). In Figs. 2-5, we observe the appearance of two disturbed plasma density minima 
under the action of ULF pulsations (the curves 1 in Figs. 2 and 3 display the action 
of the centrifugal and gravitational forces). This two minima effect could be verified by 
observations. 
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Figure captions 



Fig. 1. The dependence of the Earth's magnetic field B on the geomagnetic latitude 
(fi in the meridional plane of the dayside magnetosphere for different L in the model by 
Antonova and Shabansky (1968). 

Fig. 2. The dependence of the relative density perturbation S ior L — 6, k — 13. The 
curve 1 corresponds to gefj 7^ 0, Bio = 0, curve 2 - (^e// 7^ 0, Biq = 10~^ G, and curve 3 - 
g,ff = 0, Bio = 10-^ G. 

Fig. 3. The same as in Fig. 2 for L = 10. 

Fig. 4. The dependence of 6 distribution on the dipole {k — 0) and non-dipole {k — 13) 
geomagnetic field for L = 8 (ggfj 7^ 0, B^q — 10~^ G). 

Fig. 5. The same as in Fig. 4 for L = 10. 



